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ABSTRACT 

We present deep radio continuum observations of the cores identified as deeply embedded 
young stellar objects in the Perseus molecular cloud by the Spitzer c2d programme at a wave- 
length of 1.8 cm with the Arcminute Microkelvin Imager Large Array (AMI-LA). We detect 
72% of Class objects from this sample and 31% of Class I objects. No starless cores are 
detected. We use the flux densities measured from these data to improve constraints on the 
correlations between radio luminosity and bolometric luminosity, infra-red luminosity and, 
where measured, outflow force. We discuss the differing behaviour of these objects as a func- 
tion of protostellar class and investigate the differences in radio emission as a function of core 
mass. Two of four possible very low luminosity objects (VeLLOs) are detected at 1.8 cm. 

Key words: Radiation mechanisms:general - ISM:general - ISMxlouds - stars: formation 



1 INTRODUCTION 

Since the advent of the Spitzer space telescope, direct infrared de- 
tection of the hot dust heated by an embedded protostar has become 
the most popular method for differentiating between starless and 
protostellar cores. However, Spitzer measurements of the luminos- 
ity distribution of protostars (Dunham et al. 2008) have aggravated 
the 'luminosity problem' first articulated by Kenyon et al. (1990). 
This problem arises from the increasing number of protostellar ob- 
jects being discovered with internal luminosities Ljut < 1 'Lr-,, and 
furthermore Lj^t < 0. 1 Lrr, in the case of very low luminosity objects 
(VeLLOs). Such low luminosities violate the lower luminosity lim- 
its set by steady accretion models for young stellar objects (Evans 
et al. 2009; Shu 1977) and imply that alternative processes, such as 
non-steady accretion (Kenyon & Hartmann 1995; Young & Evans 
2005; Enoch et al. 2007) are required to explain this discrepancy. 
This luminosity problem is most difficult to rectify in very low lu- 



* We request that any reference to this paper cites "AMI Consortium: Scaife 
etal. 2011". 

t email: ascaife@cp.dias.ie 



minosity objects (VeLLOs; Young et al. 2004; Dunham et al. 2008) 
with extreme luminosities Lj^t ^0.1 Lq. The nature of these ob- 
jects is unclear, whether they are young Class protostars which 
are just powering up, or are more evolved but in a low accretion 
state (Dunham et al. 2008; Evans et al. 2009). 

Although the method is extremely reliable, Spitzer identifica- 
tion is not completely certain (Hatchell & Dunham 2009) as the 
response to extra-galactic radio galaxies is known to mimic that 
of embedded protostellar cores and can cause false detections. In 
addition, a protostellar spectrum alone will not demonstrate that a 
core is truly embedded and, like many surveys, Spitzer is also flux- 
limited, with a luminosity completeness limit of 0.004(rf/140 pc)^ 
Lp, (Dunham et al. 2008; Harvey et al. 2007). Correctly determin- 
ing the relative numbers of starless and protostellar cores in star- 
forming regions is essential for inferring timescales for the different 
stages of protostellar evolution; for low-luminosity objects it is also 
necessary in order to determine the extent of the luminosity prob- 
lem. Identifying embedded protostars by detecting their molecular 
outflows, either via high-excitation jet interactions (Herbig-Haro 
objects and shock-excited H2; eg. Davis et al. 2008; Walawen- 
der et al. 2005) or through low-excitation molecular lines such as 
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CO, can immediately identify a source as an embedded protostar, 
avoiding the contaminants of infrared colour selection. However, in 
crowded star-formation regions, confusion from neighbouring out- 
flows can frequently be an issue, as we discuss later. Outflows are 
particularly interesting for VeLLO scenarios by identifying if the 
protostar has shown more active mass ejection in the past, as mo- 
mentum deposited at earUer times remains visible in the molecular 
outflow. 

A further method for detecting protostars is through their ra- 
dio emission (see e.g. Andre, Ward-Thompson & Barsony 1993). 
This radio emission provides not only a detection method, but also 
a potential mechanism for distinguishing protostellar class via its 
correlations with other physical characteristics. Radio follow-up 
of the Dunham et al. (2008) Spitzer catalogue of low luminos- 
ity embedded objects at 16 GHz (AMI Consortium: Scaife et al. 
2010) has provided clear evidence for distinct trends in the be- 
haviour of the radio luminosity of these sources when compared 
with their bolomctric luminosity, IR luminosity and outflow force 
where molecular outflows have been measured. Similar trends are 
well-documented at lower radio frequencies for higher luminos- 
ity sources (< 10^ L©), see e.g. Anglada (1995) or Shirley et al. 

(2007) . In addition, the relationship between radio luminosity and 
the IR luminosity of these objects has shown that it may be possi- 
ble to use the radio emission as a proxy for the internal luminosity 
of low luminosity objects, a quantity which can otherwise only be 
derived through complex modelling of the IR spectra. 

Here we present new data following on from the initial work of 
AMI Consortium: Scaife et al. (2010), hereafter Paper I. The can- 
didates identified by Spitzer data in the catalogue of Dunham et al. 

(2008) were ranked as belonging to one of six "groups", with Group 
1 being those most likely to be true embedded objects, and Group 
6 those least likely. These new data cover all objects classified as 
Group 1-3 in Dunham et al. (2008) at 5 15° not covered by Pa- 
per I. These sources lie predominantly in the area of the Perseus 
molecular cloud and significantly increase the available high fre- 
quency radio data for this region. The Perseus molecular cloud is 
an exceptionally well studied star formation region with extensive 
complementary multi-frequency data available (e.g. Hatchell et al. 
2005; 2007a;b; Hatchell & Dunham 2009). The sources investi- 
gated here all form part of the low- and very low luminosity object 
population which has aggravated the 'luminosity problem'. As such 
the detailed correlations and derivations of their physical properties 
are of particular interest. 

The organization of this paper is as follows. In § 2 we describe 
the sample of targets to be observed, and in § 3 we describe the 
AMI Large Array (AMI-LA) telescope, the observations and the 
data reduction process. In § 4 we comment upon the results of the 
observations and compare them to predictions. The nature of the 
radio emission and the correlations derived in § 5 are discussed in 
§6. 



2 THE SAMPLE 

Previous centimetre-wave radio continuum follow-up (Paper I) of 
the Spitzer c2d catalogue of deeply embedded protostars (Dunham 
et al. 2008) covered all those objects which lay within clouds tar- 
geted by the AMI-SZA spinning dust sample (Scaife et al. in prep). 
Candidates for embedded objects from this catalogue were ranked 
by Dunham et al. (2008) as belonging to one of six "groups", with 
Group 1 being those most hkely to be true embedded objects, and 
Group 6 those least likely. The sample observed in this paper targets 



Table 2. AMI-LA frequency channels and primary calibrator flux densities 
measured in Jy. 
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Freq. [GHz] 


13.88 


14.63 


15.38 


16.13 


16.88 


17.63 


3C48 


1.89 


1.78 


1.68 


1.60 


1.52 


1.45 


3C286 


3.74 


3.60 


3.47 


3.35 


3.24 


3.14 


3C147 


2.85 


2.71 


2.57 


2.45 


2.35 


2.24 



all the remaining candidates not observed in Paper I and identified 
as belonging to Groups 1-3 in Dunham et al. (2008), although we 
note that no Group 2 candidates are present, above 8 = 15°. This 
declination limit is a consequence of the observing range of the 
AMI telescope; below this declination significant interference from 
geostationary satellites is experienced. The objects in this sample 
are predominantly found in the Perseus molecular cloud (22 cores) 
with a small number being located elsewhere (6 cores). We iden- 
tify the candidates from this catalogue by their catalogue number, 
i.e. [DCE08]-nnn, in column [1] of Table 1. The coordinates of the 
individual candidates are listed in columns [2] and [3] of Table 1, 
along with their Group and infrared luminosity from Dunham et al. 
(2008) in columns [4] and [5]. They have been cross-referenced 
with the sub-mm star formation survey of the Perseus molecular 
cloud (Hatchell et al. 2005; Hatchell et al. 2007a, hereafter H07), 
the designations from which are listed in column [7] and the Class 
assigned to each object from those works is listed in column [8]. 



3 OBSERVATIONS 

AMI comprises two synthesis arrays, one of ten 3.7 m antennas 
(SA) and one of eight 13 m antennas (LA), both sited at the MuUard 
Radio Astronomy Observatory at Lord's Bridge, Cambridge (AMI 
Consortium: Zwart et al. 2008). The telescope observes in the band 
13.5-17.9 GHz with eight 0.75 GHz bandwidth channels. In prac- 
tice, the two lowest frequency channels (1 & 2) are not generally 
used due to a lower response in this frequency range and interfer- 
ence from geostationary satellites. 

Observations of the 28 objects listed in Table 1 were made 
with the AMI-LA between luly 2010 and August 2010. Each tar- 
get was observed as a single pointing, with the exceptions of 
[DCE08]-048 which is located within the [DCE08]-049 pointing 
(separation's 4"), and [DCE08]-107 which is within the [DCE08]- 
108 pointing (separation's 6"). 

AMI-LA data reduction is performed using the local software 
tool REDUCE. This appUes both automatic and manual flags for in- 
terference, shadowing and hardware errors, Fourier transforms the 
lag-delay correlator data to synthesize frequency channels and per- 
forms phase and amplitude calibrations before output to disc in uv 
FITS format suitable for imaging in AIPS^ Flux (primary) calibra- 
tion is performed using short observations of 3C286 and 3C48. We 
assume I-l-Q flux densities for this source in the AMI LA channels 
consistent with the updated VLA calibration scale (Rick Perley, pri- 
vate comm.), see Table 2. Since the AMI-LA measures I-i-Q, these 
flux densities include corrections for the polarization of the calibra- 
tor sources. A correction is also made for the changing interven- 
ing air mass over the observation. From other measurements, we 
find the flux calibration is accurate to better than 5 per cent (AMI 

' http : / /www . aips . nr ao . edu/ 
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Table 1. The AMI-LA sample of embedded protostellar sources selected from the catalogue of Dunham et al. (2008). Columns are [1] source number from 
the Dunham et al. (2008) catalogue; [2] Right Ascension of source in J2000 coordinates; [3] Declination of soiffce in J2000 coordinates; [4] candidate Group 
from Dunham et al. (2008); [4] IR luminosity from Dunham et al. (2008); [5] distance in Idloparsecs; [6] cross-identification with the catalogue of HatcheU 
et al. (2007); and [7] protostellar Class, where 'S' indicates a starless core. 
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Consortium: Scaife et al. 2008; AMI Consortium: Hurley-WaUcer 
et al. 2009). Additional phase (secondary) calibration is done us- 
ing interleaved observations of calibrators selected from the Jodrell 
Bank VLA Survey (JVAS; Patnaik et al. 1992). After calibration, 

the phase is generally stable to 5° for channels 4-7, and 10° for 
channels 3 and 8. The FWHM of the primary beam of the AMI LA 
is ^ 6arcmin at 16 GHz. The data in this paper were taken with 
the AMI Large Array (AMI-LA) using channels 4 — 7, which have 
superior phase stability. Therefore measured flux densities in what 
follows are made with a total bandwidth of 3 GHz. 

Reduced data were imaged using the AIPS data package. 
CLEAN deconvolution was performed using the task IMAGR which 
applies a differential primary beam correction to the individual 
frequency channels to produce the combined frequency image. In 
what follows we use the convention: 5v °= v", where 5v is flux 
density (rather than flux, Fy = vSy), v is frequency and a is the 
spectral index. All errors quoted are 1 a. 



4 RESULTS 

The observations towards the 28 targets listed in Table 1 and de- 
scribed in § 2 are stmmiarized in Table 3. The details of each 
individual observation including both the primary and secondary 



calibration sources are listed. Maps were made using naturally 
weighted visibilities to ensure optimal noise levels, except in those 
cases where sources were separated by less than one synthesised 
beam. In these instances uniform weighting was used in order to 
improve resolution and separate the sources. The resulting rms 
noise level in each map and the dimensions of the synthesized beam 
are listed in Table 3, and where uniform weighting has been em- 
ployed it is indicated in the table. The rms noise level varies be- 
tween fields due to the different levels of data flagging required 
following periods of poor weather conditions or interference from 
non-astronomical sources, such as geostationary satellites. 

Source detection for these data was performed in the un- 
primary-beam-corrected maps, where wc identified all objects 
within the FWHM of the AMI-LA primary beam with a peak flux 
density > Sc^s, as being true sources. A full list of the sources 
detected in these fields is given in Tables Al & A2. The peak and 
integrated flux densities listed in both Table 4 and Tables Al & A2 
were extracted from the primary-beam-corrected maps. 

Integrated flux densities were determined by fitting a multi- 
variate Gaussian and a background level to each source using the 
standard AIPS task JMFIT. In addition, those objects which were 
poorly fitted by a Gaussian source model have integrated flux den- 
sities found using the fitflux program (Green 2007). This method 
calculates flux densities by removing a tilted plane fitted to the lo- 
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Table 3. AMI LA observations. Sources found within the Perseus region 
are Usted first, with those outside Perseus continued below the double line 
division. Columns are [1] soitfce number from the Dunham et al. (2008) 
catalogue; [2] AMI-LA flux calibrator; [3] AMI-LA phase calibrator; [4] 
rms noise measured from recovered map; [5] major axis of AMI-LA syn- 
thesized beam; [6] minor axis of AMI-LA synthesized beam. 
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cal background and integrating the remaining flux density. This is 
done by drawing a polygon around the source and fitting a tilted 
plane to the pixels around the edge of the polygon. Where an edge 
of the polygon crosses a region confused by another source, the 
background is subjective and this edge is omitted from the fit- 
ting. Since the extracted flux density is dependent to some de- 
gree on the background emission, we repeated this process using 
five irregular polygons, each varying slightly in shape. The final 
flux density is the average of that extracted from these five poly- 
gons. The error on these integrated flux densities is calculated as 

CJS = \J (0.05Sint)^ + ct|j + <Jj2^s, where cTfit is the standard devia- 
tion of the flux densities found from the five polygonal apertures, 
CTnns is the rms noise determined from the maps and O.OSSjnt is a 
conservative 5 per cent absolute calibration error. In the case of 
sources fitted using JMFIT Ofjt is the fitting error returned by the 
task. Where the 5 per cent calibration uncertainty is not dominant, 
errors determined for this sample are in general heavily dependent 
on the rms noise, <7nns> rather than the fitting error, cjgt, which is 
found to be small. 

The source list in Tables Al & A2 was cross referenced with 
the literature and where the position of a source coincides with a 
known object this has been indicated in Column [7]. For this identi- 
fication we employed a search radius of 30", although we note that 



most associations are within < 2". The only source significantly 
separated from its association is [DCE08J-068 and this is com- 
mented upon in later discussions. Sources are identified as extra- 
galactic radio sources where they have a counterpart in the NVSS 
(Condon et al. 1998) catalogue, and are denoted "radio". 

4.1 Expected contamination by extragalactic radio sources 

At 16 GHz we expect a certain number of extragalactic radio 
sources to be seen within each of our fields. Following Paper I, 
to quantify this number we use the 15 GHz source counts model 
from de Zotti et al. (2005) scaled to the IOC survey source counts 
(AMI Consortium: Davies et al. 2010). The average rms noise from 
our datasets is ~ 25/xJybeam~' and from this model we predict 
that we should see 0.07 sources arcmin"^, or ~ 2 radio sources 
within a 6 arcmin FWHM primary beam above a 5 a flux density 
of 125 ^Jy. Within the Perseus region we would therefore expect 
to detect 26 ± 5 extragalactic radio sources, and in the remainder 
of the sample a further 9 ± 3 sources. Tables Al & A2 shows the 
full source list for objects detected at > 5 <t in the fields. Making the 
assumption that all sources which cannot be identified with a previ- 
ously known protostellar object are extragalactic, we find 21 radio 
sources in the Perseus fields, consistent with our prediction. In the 
fields outside the Perseus region we find only 3 non-protostellar 
sources. Although this is low compared to the prediction we sug- 
gest that this is simply because the covered sky area is small and 
that this is a statistical outlier, consistent at 2 a. Of these sources, 
seven are matched to known radio sources in the NVSS catalogue 
(Condon et al. 1998) at 1.4 GHz and six also have identifications 
in the WENSS catalogue at 327 MHz (de Bruyn et al. 1999). Using 
this data we can identify all of these objects as steep spectrum radio 
sources with an average spectral index of 00327 = —0.88 ±0.15. 
Assuming a canonical spectral index for non-thermal emission of 
a = —0.7, and taking the completeness limit of the NVSS survey 
as 3.4 mjy (Condon et al. 1998) we find that all but three of the 
remaining unmatched radio sources have integrated flux densities 
below this threshold when extrapolating from 16 GHz. 

Of the three sources which do not lie below the NVSS de- 
tection threshold, one is extended and has a much lower peak flux 
density making it unlikely to have been detected in NVSS. The re- 
maining two are AMI-09 and AMI-44; both sources have flux den- 
sities above 1 mJy at 16 GHz and are relatively compact. AMI-09 is 
identified as VLA-10 in Rodriguez et al. (1999) who surveyed the 
HH 7-11 region at 3.5 cm with the VLA. We also identify it with 
[EDJ09]- 176 from Evans et al. (2009) who classify it as a YSO, and 
[GMM08]-3 from Gutermuth et al. (2008) who identify this YSO 
as an embedded Class I source. AMI-44 is also tentatively iden- 
tified as a Class I protostar in the IC 348 nebula from the Spitzer 
catalogue of Muench et al. (2007). 



5 PROTOSTELLAR SOURCES 
5.1 Perseus Region 

Those fields which lie within the Perseus region surveyed in the 
sub-irmi by H07 contain a number of known protostellar sources 
in addition to those selected from the Spitzer catalogue of DCE08. 
Table 4 lists all the protostellar sources within these fields from 
H07 and identifies them with sources detected in the AMI-LA cat- 
alogue in this work. Where a known object is not detected by the 
AMI-LA an upper limit on the 1.8 cm radio flux density of that 
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22''30™30= 30™00^ 29™30^ 29™00- 

Right Ascension (J2000) 



DCE08-043 

Figure 1. Sources lying outside the Perseus region. Greyscale and contours at 16 GHz from the AMI-LA are shown with no correction for the primary beam 
attenuation. Contours are shown at 5, 6, 7 (Jrms etc., where values for arma for each field can be found in Table 3. The position of each source from Table 5 is 
shown as an unfilled square. The IWHM of the AMI-LA primary beam is shown as a circle and the synthesised beam for each map as a filled eUipse in the 
bottom left comer. A fuU set of maps for the remainder of the sources observed in this work can be foimd in Appendix B. 
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Table 4. Physical parameters of protostellar sources within Perseus. Columns are [1] source number from the H07 catalogue; [2] source number from the 
DCE08 catalogue; [3] protostellar Class as identified in H07; [4] flux density predicted from vibrational dust emission at 16 GHz; [5] flag on presence of radio 
counterpart; [6] 1.8 cm radio flux density in excess of the thermal greybody emission measured by the AMI-LA; [7] envelope mass from H07 corrected to 
D = 250 kpc; [8] Bolometric temperature from H07; [9] Bolometric luminosity from H07 corrected to D = 250 kpc, unless otherwise indicated; [10] Outflow 
force from H07b corrected to D = 250 kpc; [11] flag on possible confusion of the outflow with other nearby objects. 
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140 


y 


499.9 


15.9 


< 25 


< 1.5 


9.6±1.0 


4 


093 





51 


y 


191.6 


9.2 


32 


1.1* 


1.4 ±0.3 


7 




I 


9 


y 


128.5 


1.7 


158 


0.8 


9.8it0.8 


10 




I 


3 


y 


578.9 


0.6 


117 


0.2 


3.4 ±0.4 


12 


105 





71 


y 


911.1 


14.0 


31 


1.8* 


4.8 ±0.8 


13 


106 





61 


y 


490.6 


6.2 


59 


1.1* 


4.0 ±0.5 


15 


104 





12 


n 


< 93 


2.7 


67 


0.32* 


0.5 ±0.2 


16/18 


107 


s 


16^ 


n 


< 110 


- 


76* 


0.18* 


- 


16/18 


108 


s 


16"^ 


n 


< 110 


- 


- 


- 


- 


27 







128 


y 


1856.8*^ 


18.8 


53 


3.3 


21.5±3.0 


28 


055 





67 


y 


1858.7*^ 


9.9 


53 


10.4 


28.4 ±2.9 


29 


056 





55 


y 


473.0 


11.1 


49 


3.8 


16.3 ±1.5 


30 







60 


y 


610.2 


10.0 


50 


2.5 


7.9 ±0.8 


31 







59 


n 


< 140 


7.3 


38 


1.2 


2.8 ± 0.5 


35 




1 


17 


n 


< 1135 


3.2 


105 


4.3 


7.3 ±0.6 


36 







31 


y 


d 


3.2 


59 


1.9 


5.1 ±0.5 


37 




I 


27 


n 


< 98 


4.3 


672 


2.0 


2.5 ±0.4 


39 


063 


I 


9(27") 


n 


< 116(98) 


1.2 


394 


1.1* 


1.4 ±0.2 


41 







231 


y 


1587.0 


31.8 


35 


5.2 


35.5 ±2.9 


42 


073 





77 


y 


391.3 


15.4 


39 


4.6 


8.1 ± 1.6 


44 







87 


y 


233.4 


16.6 


58 


13.8 


56.7 ±4.4 


48 







36 


n 


< 236 


5.0 


49 


2.0 


2.1 ± 1.6 


49 




I 


28 


n 


<217 


4.5 


140 


4.1 


H 


50 




I 


41 


n 


<245 


4.4 


88 


3.4 


129.6 ± 10.' 


52 







33 


y 


_d 


3.7 


57 


2.7 


25.4 ±2.6 


65 


071 





12(33") 


y 


227.7(206.7) 


1.2 


56 


0.24* 




69 




I 


3 


n 


< 116 


0.4 


436 


0.02 


H? 


71 


065 





12 


y 


275.6 


1.3 


37 


0.22* 


H? 


74 


064 


I 


3 


n 


<87 


0.4 


144 


0.20* 


H 


75 







4 


n 


<96 


0.7 


67 


0.06 


H 


76 


088 





104 


y 


320.2 


9.6 


109 


1.2* 


0.9 ±0.5 


77 


084 





19(37") 


y 


305.1(287.1) 


3.8 


60 


1.1* 


3.7±0.5 


80 


060 


I 


8(4*) 


n 


< 87(91) 


1.5 


86 


1.1* 




84 







3 


n 


<223 


0.5 


74 


0.04 




101 




1 


1 


y 


227.8 


0.2 


463 


1.0 


H 


Bolo62 


081 





2 


n 


<98 


0.3 


103 


0.04 


H" 




068 





4 


n 


<90 




64* 


0.45* 






080 




10+(6'') 


y 


90.1 




34* 


0.64* 






090 




3 


y 


114.0 




114* 


0.57* 






092 






y 


206.7 




47* 


0.29* 






109 




IV(9'') 


n 


< 125 




345* 


0.26* 





Note: Where no outflow force is given an "H" indicates that an outflow has been detected by harp, and "H?" indicates a possible outflow (HatcheU et al. 
2007b) but no outflow force value is available, a suffix of H" indicates the lower line-wing detection criterion (HatcheU & Dunham 2010). 
" Improved fit to modified greybody with /3 = 1 .5. 

* Improved fit to modified greybody with j3 = 2.0. 
Unresolved from neighbouring source. 

Flux density confused by neighbouring source, 
t Greybody fitted using flux densities from Di Francesco et al. (2008). 

* From DCE08. 



source is given. These limits take into account the position of the 
source within the AMI-LA primary beam. Since source detection 
is performed in the maps before correction for the primary beam 
this means that the true limit on the unattenuated flux density will 
be given by S^^ < 5 t7ttiM(''), where A(r) is the primary beam at- 
tenuation at a radial distance r from the pointing centre. 



A number of the physical properties listed in Table 4 including 
the radio luminosity depend on the distance to Perseus. Although 
many studies assume a distance for Perseus of 320 kpc (e.g. H07), 
we here assume a distance of D = 250 kpc in order to remain con- 
sistent with DCE08. Where necessary, physical quantities from the 
literature have been corrected accordingly. 
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Table 5. Physical parameters of protostellar sources outside Perseus. Columns are [1] source number from the Dunham et al. (2008) catalogue; [2] association; 
[3] protostellar Class; [4] flux density predicted from vibrational dust emission at 16 GHz; [5] flag on presence of radio counterpart; [6] 1.8 cm radio flux 
density in excess of greybody emission measured by the AMI-LA; [7] envelope mass; [8] Bolometric temperature; [9] Bolometric luminosity from DCE08; 
[9] references for physical properties and class. 



[DCE08] 


Other names 


Class 


Spred,16 


radio det. 


5l6 


Menv 


^bol 


^bol 


refs. 








(AiJy) 


(y/n) 


(lUJy) 


(Me) 


(Lo) 


(K) 




005 


IRAS04381+2540 


1 


102 


y 


637.2 


2.8 ± 0.05 


0.67 ±0.02 


131 


K08; DiF()8 


044 


L1251-AIRS3 





38 


y 


119.5 


11.09 it 0.44 


0.8 


24 


LIO; W07 


045 


L1251-AIRS4 





4 


y 


224.5 


16.46 ±0.68 


0.6 


21 


LIO; W07 


048/049 


L1251-BIRS1/2 


OH" 


141 


y 


1135.7 


17.2* 


10/0.6 


87/140 


L06; DiF08; W07 


043 


LI 251 -A 1RS2 


1 


23 


n 


< 125 


14.15±2.10 


0.4 


230 


LIO 



References: Kauffmann et al. 2008 (K()8); Di Francesco et al. 2008 (DiF08); Lee et al. 2010 (LIO); Wu et al. 2007 (W07); Lee et al. 2006 (L06). 
" as classified by this work, see text for details, 
as calculated by this work, from 850 /im flux densities in Lee et al. 2006. 



5.2 Additional targets 

There are six targets in our sample which lie outside the Perseus re- 
gion. The first of these objects, [DCE08J-005, hes within the Tau- 
rus molecular cloud and is also known as IRAS04381-h2540. We 
use the physical parameters derived for this core by Kauffmann 
et al. (2008) from MAMBO 350 ^m data. The five remaining tar- 
gets are all found within the L1251 dark cloud. [DCE08]-043, 044 
and 045 are found in the vicinity of Core A, and [DCE08]-048 and 
049 in the vicinity of Core B. Their precise identifications are listed 
in Table 5. [DCE08]-048 and [DCE08]-049 are unresolved by the 
AMI-LA and are therefore listed together. Lee et al. (2006) found 
that the these sources were both borderline Class O/I based on their 
bolometric temperatures and the ratio of the sub-nmi to bolometric 
luminosity. By also considering the ratio of 3.6 ;itm to 850 pLm flux 
density we classify [DCE08]-048, identified with L1251-B IRSl, 
as Class 0; and [DCE08J-049, identified with L1251-B IRS2, as 
Class I. The physical parameters and classifications of those ob- 
jects in L1251-A are taken from Lee et al. (2010). 



5.3 Expected vibrational dust contribution 

For each of the sources Usted in Table 4 we constrain the contri- 
bution of thermal dust emission to the measured flux density at 
16GHz by fitting to sub-mm data at 1.1 mm (Enoch et al. 2004), 
and 850 and 450 fim (Hatchell et al. 2007). For the sources in Ta- 
ble 5 we use sub-mm data from the references in that Table. Fol- 
lowing Paper I we fit modified greybody spectra of the form 

5v = V^Sv(rd), (1) 

where V is the frequency of the data point, j3 is the opacity index 
for which we assume a canonical value of j6 = 1 .8, By is the Planck 
function for a dust temperature T^. A value of j6 = 1.8 agrees with 
that used by Hatchell et al. (2007) and is close to the theoretically 
derived value of JSqhs = 1-85 found for coagulated grains with thin 
ice mantles between 350 /im and 1.3 mm (Ossenkopf & Henning 
1994). 

Tables 4 & 5 list the predicted flux density due to thermal dust 
emission at 16 GHz from these fits. Where a value of ^ = 1.8 does 
not provide a good fit to the data we repeat these fits using alternate 
values of j3 = 1.0, 1.5 & 2.0. Where appropriate a second predic- 
tion for the 16 GHz flux density due to thermal dust and based on 
the minimum value of J3 is also shown in Table 4. In practice 
a value of J3 = 1.0 is not preferred for any of our sources. This is 



perhaps unsurprising as the sample are predominantly young pro- 
tostars with low values of T^, and as such are expected to have 
larger values of jS in general (Planck Collaboration: Douspis et al. 
201 1). For those objects which are only detected at one of the three 
sub-mm wavelengths we simply extrapolate to 16 GHz using the 
canonical J3. We note that there is a degree of uncertainty in these 
predictions, largely due to the systematic uncertainties in sub-mm 
flux densities, which also in general limit the determination of J3 
(Shirley et al. 201 la). In addition, a single value of /3 may not hold 
to cm-wave frequencies, notably in the case of objects with disks 
(Shirley et al. 201 la), but we consider this model to be a reasonable 
approximation in the absence of additional sub-mm and mm-wave 
data. 

In what follows we use the measured flux density at 1.8 cm 
with the predicted greybody contribution subtracted to calculate the 
radio Ituninosity. This is to ensure that the values used are represen- 
tative only of the radio emission and do not include contributions 
from the thermal dust tail which varies greatly between sources 
(see Tables 4 and 5) and which might therefore influence any con- 
clusions being drawn from the distributions examined in the later 
stages of this paper. 

The average ratio of predicted thermal dust emission to resid- 
ual radio emission is approximately 14%. This shows that although 
the contribution is small it is not negligible. In addition we can 
see that the ratio is lower for detected Class I sources (9%) com- 
pared with Class sources (18%). This might be expected due to 
their comparatively smaller dust envelopes, but could also indicate 
an increase in radio emission relative to thermal dust emission for 
these objects. 



5.4 Luminosity and Outflow Force Correlations 

A detailed investigation of the possible causes for a correlation of 
the radio luminosity with the bolometric luminosity, infrared lu- 
minosity and outflow force was presented in Paper I and references 
therein. We do not revisit these discussions here but utilise the com- 
bined data from this work and Paper I to redraw these correlations. 
Values for the physical characteristics of these sources are listed in 
Table 4. In the case of outflow force the letter 'H' indicates that an 
outflow has been detected with harp, but that no value for its mo- 
mentum flux is available in the literature. The correlations between 
these parameters and the radio luminosity are shown in Figs. 2, 3 
and 4. 

The correlation of radio luminosity with outflow force is 
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shown in Fig. 2. Although it is evident that no improvement in this 
correlation is made by the addition of these new measurements, this 
data is consistent with the loose trend measured in Paper I. As dis- 
cussed at length in Paper I the lack of a well-defined correlation 
in this case may be due to the errors inherent in measurements of 
outflow force. 

Although Class I objects are generally considered to have 
weaker outflows than Class objects there is no clear evolution- 
ary division in the sources sampled here. Indeed the outflow force 
values for the Class I objects in this sample span almost as broad 
a range as those from Class 0. As in Paper I we plot the minimum 
outflow force required to explain the cm-wave emission as arising 
from shock ionisation (Anglada 1995); see Fig. 2. Although the 
data seems to follow the general trend described by this model not 
all objects lie above the limit indicated for an efficiency of 77 = 1, 
indicating that their radio emission cannot arise, at least solely, 
from shock ionisation. 

The correlation of radio luminosity with bolometric luminos- 
ity is shown in Fig. 3. All sources for which a bolometric luminosity 
exists in the literature are included, see Tables 4 and 5 for values 
and references. The data from this work are very consistent with 
those from Paper I. The improved correlation combining the data 
from this paper with that of Paper I, shown in Fig. 3, is 

log[ii.8cm(niJykpc2)] = -(1.74±0.18) 

+ (0.51 ±0.26) log[Lboi(L.)]. (2) 

The correlation of radio luminosity with infra-red luminosity 
is shown in Fig. 4. The data shown are those for which a source 
is found in DCE08, and which have an IR luminosity derived in 
that work. The correlation with IR luminosity found in Paper I, and 
shown as a solid line in Fig. 4, is a good fit to the combined data 
sets. Notable in this plot are the possible VeLLO sources [DCE08]- 
065 and [DCE08]-073, which both have radio luminosities in ex- 
cess of the general trend. A further object which does not follow 
the correlation is [DCE08]-055, however there is an obvious ex- 
planation for this behaviour. Although we associate [DCE08]-055 
with [H07]-028 due to proximity, the radio detection is in fact also 
completely unresolved from [H07]-027. Since there is no infra-red 
luminosity available in the literature for [H07]-027 we describe the 
LiR of [DCE081-055 ([H07]-028) as being a lower Umit, and ex- 
clude it from the fitted data shown in Fig. 4. 

The improved correlation combining the data from this paper 
with that of Paper I, is 

log[ii.8cm(mJykpc2)] = -(1.23 ±0.39) 

±(O.59±O.34)log[LiR(L0)]. (3) 



5.5 Correlation with Envelope Mass 

For each of the objects cross-identified with the 850 /im Scuba cat- 
alogue of H07 we can also investigate the correlation of the ra- 
dio luminosity with core envelope mass. These masses were deter- 
mined using the 850 ^im flux density in the following way: 

M=^^^, (4) 

making the assumptions that the dust temperature, 7^ = lOK and 
the dust opacity, K^so^m ~ 0.0012 m^ kg^' . The opacity is known 
to follow the distribution Ky = Ko(v/vo)^, where j8 = 2. 

Reversing this relationship we can see that luminosity (SyD^) 
is a linear function of core envelope mass, M. Using the envelope 
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Figure 2. Correlation of 1.8 cm radio luminosity with outflow force. Data 
from Paper I are shown as unfilled circles. Class objects from this work are 
shown as filled circles and Class I objects as stars. The dashed lines show 
the theoretical relationship between 1.8 cm radio luminosity and outflow 
force (Curiel et al. 1989) for an efficiency of T] = 1 (minimum required 
force; lower Une) and an efficiency of tj = 0.1 (upper line). 
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Figure 3. Correlation of radio luminosity with bolometric luminosity for 
objects detected at > 5 (7. Class sources from this work are shown as filled 
circles. Class I sources from this work are shown as stars, and sources from 
Paper I are shown as unfilled circles. The best fitting correlation is shown 
as a solid line. 

masses derived from the 850 /im measurements and the known be- 
haviour of KV we can predict the radio luminosity from the dust 
mass at 1.8 cm. This prediction is shown as a dashed Une in Fig. 5. 
The predicted luminosities for each object at 16 GHz based on a 
greybody fitted to sub-mm data, the flux densities for which are 
listed in Column [4] of Tables 4 and 5, are shown as unfilled squares 
in this plot and the radio luminosities in excess of these predic- 
tions, i.e. the difference between the total radio luminosity and 
the greybody prediction, are shown as filled circles and stars for 
Class and Class I objects, respectively. The sub-mm greybody 
predictions for the sample are well fitted by the mass prediction 
model, with perhaps some evidence for a lower average value of 
KV, which is known to vary, or alternatively a lower value of the 
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Figure 4. Correlation of radio luminosity with infra-red luminosity for ob- 
jects detected at > 5 c Class sources from this work are shown as filled 
circles, Class I sources from this work are shown as stars, and sources from 
Paper I are shown as unfilled circles. The best fitting correlation is shown 
as a solid line. 



emissivity index j5 used to calculate the greybody emission. For 
a detailed comparison of current measured opacities for low mass 
cores we refer the reader to Shirley et al. (201 lb). The value of 
)fg50jum used here agrees with that derived by Shirley et al. (201 lb) 
for the low mass Class core B335 who found a range of values 

K850Mm = (1.18- 1.77)+0:36 X IQ-^ m2kg-l. 

It is clear from Fig. 5 that a correlation also exists between en- 
velope mass and radio luminosity, with a Pearson correlation coef- 
ficient of r = 0.89 indicating a strong positive correlation. However, 
this correlation does not seem to hold as strongly, if at all, for the 
detected Class I objects in this sample which are all found at the 
low-mass end of the relationship and deviate from the trend. Un- 
like the luminosity and outflow force correlations this trend seems 
to show a dependency on protostellar evolution, with all Class 1 
objects lying at the low mass end. The fact that the Class I objects 
also deviate from the general correlation seen in the Class sources 
suggests, similarly to the detection statistics (see § 5.7), that the ra- 
dio emission from Class 1 objects is heavily influenced by environ- 
mental factors rather than arising as a consequence of an intrinsic 
mechanism. 

The correlation fitted to the Class sources in this sample is 



log[ii.8cm(mJykpc2 



-(2.23 ±0.65) 

+(0.68 ± 0.62) log[Menv(MQ) 



(5) 



and is plotted as a solid line in Fig. 5. 

It is known that a correlation exists between bolometric lu- 
minosity and core mass (see e.g. Planck collaboration: Montier 
et al. 2011) with the form logLboi = logA + 0.671ogMenv, where 
the constant A depends on the surface density. Combining this re- 
lationship with that fitted between radio and bolometric luminos- 
ity in this work, Eq. 2, we can predict a weak correlation be- 
tween radio luminosity and core mass with the form logLjad = 
logA' + 0.31 logMenv. The relationship shown in Equ. 5 is broadly 
consistent with this correlation due to the large scatter in the data, 
and is shown in Fig. 5(b) as a dotted line. However, it is unclear 
why this relationship should fail to apply in the case of Class I ob- 
jects. It can be seen from Fig. 3 that the bolometric luminosities 



of the Class I sources in this sample show no deviation from the 
general trend, nor with any evolutionary distinction evident in their 
distribution. 

As sources evolve through the Class and Class I stages, both 
the envelope mass and the average mass loss rate in the outflow de- 
crease (Bontemps et al. 1996; Ladd & Fuller 2002; Arce & Sargent 
2006). The observed correlation of radio luminosity, Ljad, with en- 
velope mass could also be seen as a correlation with mass loss rate. 

Models of radio emission from both ionised spherical stel- 
lar winds (Panagia & Felli 1975) and collimated winds (Reynolds 
1986) both predict a correlation between radio luminosity, L^^, and 
the rate of stellar mass loss, M. In the canonical spherical case with 
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10-5 Mgyr-l 



4/3 

(6) 



(Panagia & Felli 1975). For the case of a collimated outflow 
Reynolds (1986) also predicted the relationship 



(7) 



where the constant of proportionality included weak dependencies 
on ionisation fraction, collimation, temperature and outflow incli- 
nation. Although the best-fit regression to the Class data shown 
in Fig. 5(a) & (b) is not exactly equal to this value, it can be seen 
that these data are not inconsistent with the predicted relationship, 
plotted as a dashed line in Fig. 5(b). 

In recent models, mass accretion is suggested to be episodic, 
raising luminosities and temperatures and switching protostars 
back and forth between Class and Class I (Dunham et al. 2010). 
If the variable accretion model were correct, and radio luminosi- 
ties did correlate with mass loss rate in the outflow, then Class I 
sources should have radio luminosities higher by an order of magni- 
tude than Class sources. We find that the small number of Class I 
sources which are detected in this sample have much higher radio 
luminosities relative to the correlation seen with envelope mass for 
Class 0. However, the high number of non-detections for Class I 
sources also suggests that mass loss rate is not the only factor con- 
trolling the radio luminosity for these objects. 

We also note that the observed correlation could also be ex- 
plained by a population of dust with a similar mass but a higher 
opacity, Kgjo^m ~ 0.006 m^ kg-'. This value is very high com- 
pared with the range of values collected by Shirley et al. (2011b) 
and, even considering the degenerate contributions of mass, opac- 
ity and temperature to this relationship, it would seem unlikely that 
the emission we see at centimetre wavelengths is not in fact free- 
free emission but thermal emission from a separate population of 
dust. Alternatively the emission could also arise as a consequence 
of rotational emission from a population of very small dust grains. 
If these VSGs were mixed through the dense dust core as a func- 
tion of density then the amount of rotational emission would also 
be a function of dust mass. However, without additional supporting 
evidence we do not explore this possibility further. 



5.6 Correlation with Bolometric Temperature 

The bolometric temperature, values taken from H07, is representa- 
tive of the dust temperature in the cores - rather than the tempera- 
ture of any ionized gas. Consequently the lack of correlation with 
radio luminosity is perhaps not surprising. 
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Figure 5. Correlation of radio luminosity with envelope mass, (a) Linear 
correlation, and (b) Logarithmic correlation. In (a) two data points are pre- 
sented for each object: unfilled circles show the sub-mm greybody predic- 
tions for 16 GHz luminosities from thermal dust emission, filled circles and 
stars show the measured radio luminosities in excess of the greybody pre- 
dictions for Class and Class I objects, respectively; a 16 GHz luminosity 
prediction from envelope mass is shown as a dashed line; in (b) the radio 
luminosity data only is plotted; a best fitting regression is shown as a solid 
line, the relationship L„(] o= M^/' is shown scaled to fit the Class data as 
a dashed hne, and the predicted relationship based on the that of radio and 
bolometric limiinosity is shown as a dotted Une, see text for details. 



On average the bolometric temperature is thought to increase 
with bolometric luminosity, however within individual classes of 
protostellar object this trend is extremely weak and indeed it is pos- 
sible that it is itself simply an effect of observational biases in the 
sub-mm (Kauffmann et al. 2008). Although MAMBO follow-up of 
the Spitzer c2d catalogue showed a correlation between bolometric 
temperature and envelope mass (Kauffmaim et al. 2008), this was 
only really evident when combining objects in Classes O, I and 11. 
Individual classes did not portray a strong trend between these pa- 
rameters and the same is true for the sample of objects considered 
here. When Classes O and I are considered together a weak nega- 
tive correlation is foimd, r = —0.41; however if the Class objects 



alone are considered then this correlation reduces even further, and 
in fact reverses, to r = 0.24. 

When comparing bolometric temperatures determined using 
IRAC data (H07) and those found from fluxes calculated using 
Spitzer SEDs (DCE08) there is a marked difference for this sam- 
ple, with those derived from IRAS being, on average, a factor of 1.6 
higher than the corresponding Spitzer values. Bolometric temper- 
ature is one of the indicators for protostellar class and within this 
sample three sources with 70 ^ Tboi ^ 150, will have this indicator 
change from Class I to when the Spitzer temperature is used. For 
two of these sources, [DCE08]-060 ([H07]-80) and [DCE08]-064 
([H07]-074), their overall classification will change from Class I 
to Class when Spitzer temperatures are utilised rather than those 
from IRAC. However, neither is detected at 1.8 cm by the AMI-LA. 
The third, [DCE08]-088 ([H07]-076), is already identified as Class 
and is detected. 



5.7 Detections and Non-detections 

From the 28 objects targeted in the Spitzer catalogue of embedded 
sources (DCE08) there are nine in the Perseus region which we do 
not detect. Of these, four are designated Group 1 by DCE08 and are 
therefore most likely to be true embedded protostars. From these 
Group 1 objects two are Class I following the revised classifications 
from § 5.6 ([DCE08]-063 & 064) and two are Class ([DCE08]- 
068 & 081). [DCE081-068 is identified as Class following an 
identification with [H07]-75 in Hatchell & Dunham (2010), how- 
ever we note that the positions of these sources are in fact offset 
by 45" making their correspondence uncertain. This may be an ex- 
ample of a false protostellar identification due to confusion with 
an external galaxy in DCE08. [DCE08]-081 was initially identi- 
fied as Class by H07, but no outflow was detected by Hatchell 
et al. (2007b). An outflow was identified by Hatchell & Dunham 
(2010) using the "lower line-wing criteria", with the caveat that it 
may have been confused with the outflow from [H07]-081. Five of 
the six Group 3 objects are not detected. Of these two are consid- 
ered starless ([DCE08]-107 & 108), one is Class I ([DCE08]-060), 
one has no classification and no counterpart in H07 ([DCE08]-109). 
The remaining non-detection is [DCE08]-104, which is a Class 
source. 

Of the sources outside the Perseus region all Group 1 objects 
are detected, although we note that [DCE08]-048 and 049 are un- 
resolved and it is likely that the radio detection is dominated by 
[DCE081-048. [DCE081-043 (Class I), flie only Group 3 object out- 
side Perseus is not detected. From these numbers the radio detec- 
tion rates are 81 % and 29 % for the Group 1 and Group 3 objects, 
respectively. 

The detection statistics for the whole sample divided by proto- 
stellar Class and candidate Group are summarised in Table 6. The 
results are very similar to those found in Paper I for VeLLO and 
Class sources. For Class I sources the detection rate is lower 
in this work compared with that of Paper I (71 %). The consis- 
tency in detection for Class sources, and lack of consistency 
for Class I sources, may indicate that the radio emission mecha- 
nisms for the two types of protostellar objects differ. Consistency 
indicates that the emission from Class objects is dependent on 
an intrinsic mechanism, whereas for Class I objects environmen- 
tal factors are important leading to different detection statistics in 
different regions. We note that a comparison of the statistics for 
VeLLO sources is not conclusive due to the small number of these 
objects in this work. We also tabulate the detection statistics for 
the much larger number of protostellar sources which lie within 



AMI-LA observations of Perseus 1 1 



our fields identified from the catalogue of H07, however we note 
that the primary beam attenuation away from the pointing centre 
will bias these detections non-uniformly. We include starless cores 
identified in H07 in these statistics, but do not list the positions for 
these objects in Table 4 as they arc uniformly not detected. In spite 
of the differences between these samples the detection statistics for 
Class and Class I objects are very similar for the original DCE08 
sample and the extended dataset. We note that these numbers in- 
clude [DCE08]-049 as a detection, although it is not possible to say 
conclusively whether it is contributing to the flux density of AMI- 
49, which combines [DCE08]-048 and [DCE08]-049. The effect of 
considering it a non-detection would be to reduce the Class I de- 
tection rate to 17(25)% for the original sample and 18(22)% for the 
extended sample. 

Since the correlation of radio luminosity with bolometric lu- 
minosity is the most well defined we examine the distribution of 
these non-detections with respect to this relationship. There are 
two distinct regions of non-detection. The first comprises low- and 
very low bolometric luminosity objects. The non-detection limits 
for these sources lie entirely above the general correlation and we 
can therefore assume that the fact they are not detected is biased by 
their low intrinsic luminosities. The second group is found at higher 
bolometric luminosities (Lbol > 1 Lr:,,) and lies predominantly be- 
low the general trend. For these objects there are two possibili- 
ties, firstly that their non-detection is simply an observational bias 
caused by the sensitivity limit of the observations, or secondly that 
they do not follow the general trend and have abnormally low, or 
zero, radio emission. There is no evident common characteristic be- 
tween these sources which might be used to distinguish them from 
the general population, and in the absence of more sensitive obser- 
vations it is not possible to say definitively if this is indeed the case 
as the sensitivity limit is not significantly removed from the general 
trend. 



5.8 VeLLOs 

There are three sources in this sample which can be defined as VeL- 
LOs based on their modelled internal luminosities by DCE08, they 
are [DCE08]-064, 065 and 081. The source [DCE08]-073 has an 
infra-red luminosity of Lir = 0.018 L,-:,, well below the rough cut- 
off of LjR = 0.05 L0 that generally indicates VeLLO status, which 
may also make it a VeLLO candidate or borderline case. Of these 
sources, we detect [DCE08]-065 and [DCE08]-073 at 1.8 cm with 
the AMI-LA. Both of the detected sources are identified as Class 
from H07. Of the undetected sources, [DCE08]-081 is identified as 
Class and [DCE08]-064 as Class I by H07, although we note that 
[DCE08]-064 would be re-classified as Class using its bolomet- 
ric temperature from DCE08 as discussed in Section 5.6. [DCE08]- 
073 and [DCE08]-065 deviate from the general correlation trend of 
radio luminosity to infra-red luminosity, having a radio luminosi- 
ties in excess of that expected, and this has already been noted in 
Section 5.4. However, with two sources of similar infrared luminos- 
ity from Paper I having radio luminosities below the general trend 
it is likely that this is simply due to scatter in an under-sampled lu- 
minosity range rather than representing a characteristic deviation. 
[DCE08]-073 is the only VeLLO object with a measured outflow 
force and we note from Fig. 2 that the radio luminosity for this 
source is consistent with having arisen as a consequence of shock 
ionisation, unlike the VeLLO source L1014 (Shirley et al. 2006) 
which has an outflow momentimi significantly too low to explain 
its radio emission. 



Table 6. Summary of detection statistics for Table 4. The original sample 
includes only objects from DCE08, listed in Table 1. The extended sample 
includes all sources identified in H07 and listed in Table 4, inclusive of the 
original sample. 



Class 


Present 


Detected 
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Original 


sample: 






VeLLO 
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50 





16(18) 


13(13) 


81(72) 


I 


6(4) 


2(2) 


33(50)t 


Starless 


2 








Group 1 


21 


17 


81 


Group 3 


7 


2 


29 


Extended sample: 









25(27) 


18(18) 


72(67) 


I 


11(9) 


3(3) 


27(33)t 


Starless 


15 









Note: Figures in brackets indicate revised numbers following 
re-classification of sources based on bolometric temperatures 
derived from Spitzer data in DCE08. 

+ By considering [DCE08]-049 as a non-detection these figures 
are reduced to 17(25)% for the original sample and 18(22)% for 
the extended sample, see text for details. 

6 DISCUSSION AND CONCLUSIONS 

The evidence for hot gas in low mass protostars that could give 
rise to radio frcc-frcc emission has been strengthened by recent 
analyses of Herschel-PACS data (van Kempen et al. 2010) which 
have revealed spectral features indicative of small-scale shocks cre- 
ated along the cavity walls of protostellar outflows. Of these, C- 
type shocks give rise to H2O emission, and J~type shocks within 
the lower density jet give rise to observed Ol features. Combina- 
tions of these shocks can also give rise to observed OH emission. 
However, the gas in such shocks is heated to temperatures of up 
to only 3000 K (van Kempen et al. 2010), just sufficient to pro- 
duce free-free emission. Typical temperatures for free-free emis- 
sion are in the range 4000 ^ Te < 20000 K, with a mean value of 
approximately 8000 K (Dickinson et al. 2003). Below 5000 K emis- 
sion through molecular lines and bands becomes significant. Indeed 
the theoretical lower limit that we show in Fig. 2 from Anglada 
(1996) for the amount of free-free emission from shock ionisation 
for a given outflow force, F, assumes Je = 10* K. This relation- 
ship can be re-expressed to include a temperature dependence as 
Lrad °= ^'outTi implying that for a lower temperature such as 
that found in van Kempen et al. (2010) a smaller amount of radio 
emission would be expected for a given outflow force. This would 
make it more difficult to explain the observed radio luminosities as 
arising from shock ionisation. 

It is unclear why free-free emission produced by shock ioni- 
sation in such regions should be correlated strongly with envelope 
mass but less so with outflow force, as has been observed in this 
work. Since the stellar mass, M*, determines the outflow velocity it 
could be supposed that emission from shock ionisation in a molec- 
ular outflow might be correlated with . If, as demonstrated by 
Chabrier & Hennebelle (2010), the stellar mass is a linear function 
of the core mass then this may lead to a correlation between en- 
velope mass and radio luminosity from shock ionisation. However, 
the formation of shocks along the length of an outflow is depen- 
dent not only on the velocity of that outflow, but also on the lo- 
cal enviroimient. To establish a clear correlation would require not 
only uniformity of the local environment but also of the nimiber 
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of shocks formed by an outflow of particular velocity. In addition 
we note that the measured outflow force for protostellar objects is 
not necessarily a current value but instead represents the averaged 
historical accretion history of the source and therefore may not nec- 
essarily correspond to the radio luminosity, which depends on the 
"current" protostellar state. The radio emission from protostars is 
also known to possess a degree of variabiUty, as observed in the 
VeLLO source L1014 (Shirley et al. 2007), which may contribute 
to the spread in the correlations presented here. The data here are 
observed over only a single epoch, and a degree of decorrelation 
due to variability will be common to all the measured correlations. 

Since the free-free emission (radio luminosity) observed here 
demonstrates a correlation with envelope mass, which is further- 
more consistent with that expected from free-free emission as a 
consequence of partially ionised stellar winds (both spherical and 
coUimated), it is possible that the thermal and non-thermal radio 
flux density (see e.g. Rodriguez et al. 1999) detected from low mass 
protostars arises as a consequence not only of different emission 
mechanisms, but also of different underlying astrophysical pro- 
cesses. In order to distinguish these mechanisms, which occur at 
separated locations within a protostellar system, much higher reso- 
lution observations are required than can be provided here. 

In addition to frec-frce and synchrotron emission, we must 
also consider the possibility that the correlation of radio luminos- 
ity with mass is a consequence of the cm-wave emission being due 
to a population of very cold dust. To check this hypothesis, the 
emission observed is fitted by a modified blackbody representing a 
possible very cold dust component. However, in order to accommo- 
date the amount of emission seen at 1.8 cm a dust temperature of 
< 5 K and a value of j3 << I are required. Such a low cold ther- 
mal equilibrium temperature for big dust grains is unrealistic, as is 
the spectral dust emissivity index. Current observations suggest /3 
to be in the range 1-2.5 (e.g. Boulanger et al. 1996; Paradis et al. 
2009) and the Kramers-Konig equations also suggest a value of 1 
as a theoretical lower limit to j8 (Emerson 1988). The existence of 
very cold dust that would explain the observed centimetre emission 
therefore seems very unlikely. 

In sunmiary we find five main conclusions from this work: 

• We find that 72% of Class objects in the extended sample 
have detected radio counterparts at 16 GHz, compared with a lower 
detection rate of 27% for Class I objects. No starless cores were 
detected. This detection rate for Class I sources differs from that 
of Paper I and we hypothesise that the discrepancy may be due 
to environmental effects causing localised differences in the radio 
emission from Class I objects, whereas radio emission from Class 
objects arises as a consequence of an intrinsic emission mechanism. 

• These new data strengthen the correlations from Paper I of 
radio limrinosity with bolometric and infra-red limiinosity in the 
low luminosity limit and significantly increase the available radio 
data for low and very low luminosity protostars. 

• We find no improvement in the correlation of radio luminos- 
ity with outflow force from these new data. We suggest that for the 
low temperatures found in shocks along molecular outflow cavi- 
ties the measured outflow forces cannot explain the observed radio 
luminosities. 

• We find a correlation of radio luminosity with core mass for 
Class objects which does not hold for Class I objects. We suggest 
that this difference is broadly consistent with theories of episodic 
mass accretion, but requires additional data for Class I sources in 
order to be confirmed. 

• We hypothesize that the observed correlation with core mass 



seen for Class objects may be due to the predicted relationship 
between stellar mass loss rate and radio luminosity, as suggested 
by Panagia & FeUi (1975) and Reynolds (1986). With the currently 
available data it is not possible to distinguish the case of an ionised 
spherical stellar wind from a collimated outflow, however this cor- 
relation suggests that such a mechanism is more likely to account 
for the observed free-free radio emission than the alternatively hy- 
pothesised mechanism of shock ionisation. 
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0, Class 1 and starless cores from Hatchell et al. (2007) are indicated 
as crosses ("+"), stars and unfilled circles, respectively. 



APPENDIX A: SOURCE LIST 

Sources detected in the AMI-LA fields. Column [1] indicates the 
AMI-LA source designation in increasing Right Ascension, col- 
umn [2] indicates the field in which the source was detected, col- 
vmin [3] lists the Right Ascension in J2000 coordinates, colunm 
[4] lists the source Declination in J2000 coordinates, column [5] 
lists the peak flux density of the source at 16 GHz in units of 
/xJybeam^' and column [6] lists the integrated flux density in 
units of mjy; column [7] lists known associations for each source. 
Sources are designated as "radio" where they have a counterpart in 
the NVSS catalogue (Condon et al. 1998). 



APPENDIX B: AMI-LA MAPS 

Maps of each field are shown. AMI-LA data is shown uncorrected 
for the primary beam response as greyscale and contours of 5, 6, 7, 
8 a etc to 1 njybeam"^ and then in increments of 1 mjybeam"^ 
The AMI-LA primary beam response is shown as a solid circle and 
the PSR for each dataset as a filled ellipse in the bottom left cor- 
ner. The positions of Spitzer embedded protostellar cores (Dunham 
et al. 2008) are shown as unfilled squares, and the positions of Class 
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Table Al. Detected sources in the Perseus region at 1 6 GHz witti ttie AMI-LA. 
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Figure Bl. 
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Table A2. Detected sources outside the Perseus region at 1 6 GHz witii ttie AMI-LA. 
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